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Abstract 

The properties of vortices in Josephson junction arrays are investigated in 
the quantum regime near the superconductor-insulator transition. We de- 
rive and study an effective action for vortex dynamics that is valid in the 
region where the charging energy is comparable to the Josephson coupling 
energy. In the superconducting phase the onset of quantum effects reduces 
the vortex mass and depinning current. In the case of long range Coulomb in- 
teraction between Cooper pairs we find that as the transition is approached, 
the velocity window in which ballistic vortex motion is possible grows. At 
the superconductor-insulator transition the vortex mass vanishes and vortices 
and spinwaves decouple. In the case of on-site Coulomb repulsion (which is 
of relevance for superconducting granular films) the vortex mass it is sample- 
size dependent in the superconducting phase, but stays finite at the critical 
point where it is scale invariant. The relation of our work to experiment is 
discussed. 

PACS numbers: 74.50, 74.60.Ge, 74.65+n 
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I. INTRODUCTION 



In classical Josephson junction arrays (JJA) vortices are well known topological ex- 
citations that characterize both the dynamical and thermodynamical properties of these 
systems. At low temperatures these excitations are bound in dipoles of opposite vorticity 
and at a critical temperature they unbind in a Kosterlitz-Thouless-Berezinskii (KTB) phase 
transition [|TJ leading the system from the superconducting (SC) to a resistive phase. 

When an external driving current is applied to the array, free vortices (as induced by 
applying a magnetic field) move and a voltage drop accross the sample appears. In classical 
JJA's, in which the Josephson coupling Ej is much larger than the charging energy Ec, 
vortex motion is diffusive. Some time ago it was realized that charging effects yield dynamics 
and a mass for the vortices in the SC phase @-fll- The electrostatic energy stored in the 
junction capacitances can be interpreted as a kinetic energy due to the vortex motion and as 
a consequence a mass can be attributed to the vortex. The mass in this regime was found to 
be proportional to the junction capacitance C and in case of small damping ballistic vortex 
motion was predicted. Modern lithographic techniques allow for the realization of JJA's 
with junction resistance much larger than h/e 2 , thereby making possible an experimental 
check of these predictions. Very recently, in an important experiment, van der Zant et al. 
H demonstrated the existence of such ballistic vortices in triangular arrays. 

The possibility of ballistic motion depends crucially on the dissipation mechanisms avail- 
able to a vortex. Nowadays it is possible to fabricate high quality JJA's in which no ohmic 
dissipation is present f|. If the vortex velocities are below a certain threshold (related to 
the superconducting gap), also the quasiparticles are frozen and do not contribute to the 
damping. It was pointed out in recent analytical [[f]J§ and numerical |J studies that another 
mechanism by which a classical moving vortex can loose its kinetic energy is by emitting 
spinwaves. Depending on the spinwave spectrum, it has been shown that in triangular arrays 
there exists a small window of velocities for the vortex to move over the pinning potential 
provided by the lattice without suffering too much damping. For square classical arrays 
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ballistic vortex motion is not possible. 

All the previous investigations were mainly concerned with the case in which the Joseph- 
son coupling Ej is much larger than the charging energy Ec (and larger than the temperature 
T). When these two energies become comparable quantum effects play a major role and, 
eventually, at a certain critical value of their ratio, a superconductor-insulator (S-I) transi- 



tion at zero temperature takes place |10Hl2|. This phase transition separates regions where 
either the Cooper pairs or the vortices are localized. In the insulating phase vortices are 
delocalized, as are the Cooper pairs in the SC phase. Thus, close to the S-I transition one 
expects vortex properties to reflect the presence and nature of the transition. One may also 
expect that experiments on vortices in this region [13] yield more insight into the phase 
transition. 

The analysis of vortex dynamics as one approaches the S-I transition will be the aim of 
this paper. We will consider a square JJA at T = with a superimposed external magnetic 
field such that a small amount of vortices of one sign are induced in the system. If the 
magnetic frustration is very small one may suppose that the vortex motion is not influenced 
by their mutual interaction. By applying an external current the dynamics can be studied. 
When Ej ^> Ec the description has been formulated in |2|-^|,|7|,|8|] . Vortex motion is then 
described by a phenomenological equation of motion for the vortex coordinate x of the form 



m 



M v 'x + r\x = 2nEjI /I cr + nUb ar sin(27rx) , (1) 

where M v is the vortex mass, rj a phenomenological damping, I the applied current, I cr the 
junction critical current, and Ub ar the height of the potential energy barrier for a vortex 
to go from one plaquette to a neighboring one in the potential landscape provided by the 
array. When we enter the quantum regime fluctuations of vortex anti- vortex pairs (dipoles) 
start to be relevant and they will interact with the moving vortex. As a result the vortex 
dynamics will be modified. 

In order to describe this regime, we first derive an effective action for one single vortex 



in the next section. We will show that it depends on the charge-charge correlation function; 
Eq. (|5[) is the central result of our work. It implies that a description in terms of an equation 
of motion is still possible, however, with modified coefficients. General expressions for the 
vortex mass and the spin wave damping that are not constrained by the inequality Ej ^> Eq 
and are also valid when the Josephson energy is comparable with the charging energy are 
then presented. After that we briefly show in Sect.|T| how to recover the previously known 
results in the limit of large Josephson coupling. Analytic approaches in the two opposite 
limits in which either the junction capacitance or the capacitance to the ground is dominant 
are discussed in Sects.|IV| and |V| respectively. For long range Coulomb interactions we find 
that closer to the S-I transition the vortex mass clearly deviates from the classical results 
and is not simply proportional to the junction capacitance. Furthermore a velocity window 
for ballistic motion opens if one approaches the transition, also for square arrays in which 
ballistic motion in the limit Ej ^> Eq is impossible. The effect of the capacitance to the 
ground has until now been overlooked for JJA's, but it gives a contribution to the vortex 
mass which is of the same order of magnitude (for a typical experimental setup) as the mass 
due to the junction capacitance and therefore, we argue, it should be taken into account. 
In the limit of short range Coulomb interactions we notice that the mass should be finite 
at the transition and the vortex dynamics is governed by the same correlation functions 
that lead to the universal conductance at the S-I transition. Therefore we feel that the 
understanding of vortex dynamics in this regime may shed light on the source of dissipation 



which is responsible for the metallic behavior at the S-I transition [|15|,[nj. Sect.|V| discusses 
Monte Carlo simulations that verify our conclusions. The main results of our work are then 
summarized in the last section. 

II. THE EFFECTIVE ACTION 

We start from the well known Hamiltonian for a Josephson junction array 
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H = \ E QiC^Qj -EjY, cos(0, - 0,) + £ k ■ V0, , (2) 

ij <«j> » 

where Qi and 0j denote respectively the charge and the phase of the superconducting order 
parameter of the i-th island. The external current is denoted by I. The typical energy scales 
are the Josephson coupling Ej and the charging energy Ec = e 2 /2C or E a = e 2 /2C a . C Q 
and C are, respectively, the ground capacitance and the nearest neighbor capacitance that 
constitute the capacitance matrix Cy. The range of the electrostatic interaction between 
Cooper pairs, which is described by the inverse capacitance matrix C -1 , is A -1 = ^JC/C Q . 
The quantum mechanical description is completed by the commutation relation [Qi, <f>j] = 
2eSiji. 



By means of well established duality transformations 17] it is possible to recast the 



partition function in terms of the topological excitations of the system. For the Hamiltonian 



under consideration this was done in Refs. |12]JI5|, where a detailed derivation is presented. 



In order to make this paper self-contained we sketch the main steps in appendix 0. The 
result is a description in terms of two discrete gasses q and v that denote charges and vortices 
respectively. In this formulation the partition function for a JJA is obtained summing over 



all configurations of charges and vortices |12 



2=EE exp{-%,v]} 

where the action is 

S[q, v] = errEj ^ v^G^Vj^ + ^^£- ^ 7,. ; /',//;.; 

ij,r 71 ij,r 

+ i ]T qi,T®ijVj, T + ie Yl h-r ' ^®ij v j,T + ^ E ]C %rGijqj, T (3) 

ijjT ijjT J ij,r 

where the kernels 0^ = arctan (f^zf 1 ) an d Gij = — In | r^ — rj | were introduced. is related 
to the capacitance matrix as Uij = 2irCCn , it has the explicit form Uy = Kq{\ | r, — r,j |), 
where Kq is a modified Bessel function [T^. Spinwaves are described by the last term in 



Eq.([3]). We introduced N T time-slices to describe the quantum dynamics of the system. The 



lattice spacing in the time direction is denoted by e and N T e = /3 is the inverse temperature 
(see App.|A]). The phase diagram for the model described by (^j was studied in ref. (12 



and exhibits a zero temperature S-I transition (an analysis based on a Kosterlitz type RG 



procedure is in progress [jnj ) . 

Eq.(£J) will be our starting point. We seek an effective action for a single vortex (with 
coordinate r(r)) that includes the effect of the interaction with fluctuating charges and 
other vortices (present in the system because of quantum fluctuations). It is obtained 
by a summation in the partition function over all configurations of the charges and other 
vortices. It turns out to be more transparent to introduce the trajectory of the vortex 
Vi tT = v5(ri — r(r)). Here v = ±1 for a vortex or anti- vortex respectively. Formally the 
effective action can now be written as 

S ef f = - In (e2irEjv^2v itT Gij8(rj -r(r)) + 
^EW^-rMj+weEi- Ve^(r,--r(r))) , (4) 

ij,r ij,r 

where the average is to be taken with the action (|3[). The first term describes the static 
interaction with other vortices, whereas the second describes the dynamical interaction with 
charges. This expression is formally exact, but difficult to evaluate because of the nonlin- 
earity of the action (|3|). As we are interested in the kinetic contribution (quadratic in the 
vortex velocities) it is, however, sufficient to expand the average in (H) in cumulants and stop 
at the second order in the vortex velocities r(r). For a uniform external current distribution 
we find 

Seff = lY,r a (T)M ab (r(r) - r(r'), r - r r )r\r') + 2mve £ e ab I a T r\r) , 

M ab = Y J V a Q(r(r) - r i )(g ir9fer ,)V 6 e(r fe - r(r')) , (5) 

jk 

where a,b = x,y and e ab is the anti-symmetric tensor. Thus, vortex dynamics is governed by 
the charge-charge correlation, which depends on the full coupled charge vortex gas (CCVG) 
Eq.(^|). The effective action Eq.fl5|) describes the dynamical vortex properties for all values 
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of Ej/Ec and is therefore a good starting point for the investigation of vortex properties 
near the S-I transition. 



terms in the cumulant expansion. This is certainly correct in the Ej ^> Eq limit where (as 
discussed in detail in the next section) the charges can be considered as continuous variables 
and vortex fluctuations can be disregarded. In general, however, the average defined in 
Eq.(f|) is far from Gaussian (q and v are integer valued fields). As long as we are in the 
superconducting phase, however, the charges are strongly fluctuating and the vortices are 
still bound in dipoles. Therefore we expect that the higher order cumulants are still not 
important A full description of the vortex motion in the resistive region, nevertheless, 
may require the analysis of a dynamical equation that contains also terms proportional to 
higher powers of the velocity. 



The expression (||) yields the known results ||[7|,|| i n the classical limit where Ej 3> Ec- 
In this region of the phase diagram the Josephson coupling dominates over the electrostatic 
energy. Far from the transition vortex fluctuations due to quantum effects are suppressed and 
they may be neglected. In this regime the charges are wildly fluctuating (as to constitute the 
supercurrents that keep the SC phases well defined) and may be considered to be continuous 
variables. 

Thus in the classical limit we may concentrate on the charge part of the action (|3|). It is 
rewritten conveniently as 



where the plasma frequency uj p = ^8EjEq was introduced. The charge-charge correlation 
is half of the inverse of the kernel Q. It is 



In order to obtain the vortex effective action in Eq.(^|) we disregarded all higher order 



III. THE CLASSICAL LIMIT 




(??)*,. 



(Ejk 2 /e)/(ul + ul), ul = ulU k /G k , 



(7) 
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The spinwave dispersion is described by uj^. It is optical, i.e. u>k = co p , for long 
range Coulomb interactions, whereas for on-site interactions we have = oJ p k. Here 



Up = y/8EjEo is the plasma frequency for the case of on-site Coulomb interactions. 

The action (|5]) reduces to that of a free particle in the limit of small velocities r(r). Since 
the charge-charge correlation (|7|) is short range in time we may put r(r) = r(r') in Eq.(^). 
The corresponding adiabatic vortex mass M v is 

M v = eJ2M xx (0,r) , (8) 

T 

which reduces in the classical limit to M v = Mes + m (-^) M a , where Mes — tt 2 /AEc 
is the Eckern-Schmid mass and M a = tc/8E . Thus both self- and nearest neighbour 
capacitances yield a contribution to the mass. The self capacitance contribution depends 
on the system size L. For generic sample sizes and capacitance ratio's the new contribution 
(which has been overlooked so far for JJA's since it makes no sense in the thermodynamic 



limit, see ref. []2T| for size dependent vortex masses in different systems) is somewhat smaller 
than the Eckern-Schmid mass. 

The instanton action Si nst , related to a hop from one plaquette to a neighbouring one, 
determines tunnel rates and the depinning current. As vortex trajectory we now take Vi T = 
Vi ;T+e — v i)T = 5 Ti t[5i ;X+a — 5i >x ] for a hop from x, t — > x + a, t + e. This may be inserted in the 
CCVG action (|3]) to find the general result 

S ins t = \M xx (0,0). (9) 

In the classical limit we recover all known results \\2,'J\\ , i.e. for general capacitance matrix 

ttE 



5, 



inst A 



(10) 



which reduces to Si nst = tt 3 ^ 2 Ej/Auj p and Si nst = tt 2 Ej/2lu p for C = and C a = re- 
spectively. The general form of an instanton action in the WKB approximation, which is 
proportional to the square root of mass times barrier height, Si nst ~ ^M v Ubar, determines 
the barrier Ub ar for a vortex to hop. The depinning current Id ep is half the barrier height. 
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Since deep in the classical limit Ubar = Ej/5 [Q, we may establish Si nst = Tc^j5M v Ubar/2 
for the C = case. 

The spinwave damping that a moving vortex experiences may also be calculated from 
(|). Varying the vortex coordinate r a (r) in Eq.(|) yields the equation of motion 

2me ab I b /I cr = U T Y,M a Mr)-r{r'),T-r')r b {r') (11) 

6 t' 

and its constant velocity solutions in the presence of an external current determine the 
nonlinear relation between driving current and vortex velocity, once the charge-charge cor- 
relation is analytically continued (i.e. sending iu u — > uj + id) to real frequencies ||. The 
relevant information is in the real part of Eq . (|TTD , which reads in Fourier components and 
for a constant vortex velocity r(r) = (v, 0) 

I y /Icr = lJd^ 1^ d 2 k^[S(u - Wfc) + S{u + u k )\8{u - vk x ) . (12) 

The delta functions express the spinwave dispersion (from the analytic continuation of the 
charge-charge correlation) and the vortex dispersion respectively. The overlap integral de- 
termines the amount of dissipation a moving vortex suffers from coupling to spinwaves. If 
we adopt the smooth momentum integration cut-off / d 2 k — > 2ti J °° dkk exp(— k/^/2n) that 
was introduced in Ref. 0, we recover in the classical limit the results of Refs. 0H, see Fig.l 
for the current vs. velocity relation for long range Coulomb interactions. Note that the min- 
imum velocity that a vortex needs to move over the pinning potential of the lattice follows 
from the phenomenological equation of motion Eq.fll]) by demanding that ^M v v 2 > Ubar and 
is about 0.14o;p (see the dotted line in Fig.l) which is also the velocity at which the spinwave 
damping sets in. Therefore ballistic vortex motion is almost impossible in classical square 
arrays. In triangular arrays, however, a similar analysis yields a somewhat wider velocity 
window || . In the next section we show how the inclusion of quantum effects contributes to 
the opening of a more robust velocity interval, also for square arrays, where ballistic motion 
can be observed. 
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IV. LONG RANGE COULOMB INTERACTIONS 



When the ratio Ej/Ec decreases the charge-charge correlation must be calculated beyond 
the classical approximation. We first consider the case of long range Coulomb interactions 
between Cooper pairs. According to the arguments given in Refs. ]II2|,[23, the zero tem- 
perature S-I phase transition is presumably of the KTB type, as is the finite temperature 
transition to the resistive phase. This means that no dimensional cross-over takes place at 
zero temperature, or in other words the dynamical critical exponent z equals zero. 

Thus we are led to conclude that the vortex fugacity in the superconducting phase scales 
to zero in the renormalization group sense also at zero temperature. This means that the 
charge-charge correlation function in the SC phase may still be evaluated in the absence of 
vortex fluctuations. Therefore we consider again only the charge part Eq.(^) of the action 
Eq.(^|), but in contrast to the classical limit we now treat the charges as discrete variables. 

The charge-charge correlation function may be rewritten as (see Appendix [B 1| for the 
derivation) the classical result minus a correction 

(q JT qkr') = ^Qjkrr' ~ Qjmrti^rntM Qnkt'r'i ( 13 ) 

— mntt' 

where the correlation function of the dual variables <ft is now to be calculated using the 



sine- Gordon-like action |25 



S[<P] = vr 2 £ £ faQTwh* ~ H £ cos(27T0 it ) (14) 

ij tt' it 

We will calculate the charge-charge correlation function in a self consistent harmonic approx- 
imation which is valid not too close to the transition point. It amounts to the replacement of 
the nonlinear cosine term in the Hamiltonian by a mass term (or inverse correlation lenght) 
which is determined selfconsistently by means of the Bogoliubov variational principle ||26|| . 
The constant H is related to the fugacity for charges in the original model, it is H" 1 = 2eEc- 
We take a trial action with the cos(27T0 i4 ) replaced by a mass term 

s[4>] = \ £ £ <t>M« 2 Qw> + M«'r%. (is) 
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and determine the correlation length £ (or inverse mass) from the selfconsistency equation 

r 2 



Ati 2 H 



exp(-27r 2 (0 2 ) ?2 ) (16) 



in the usual way 28]. The result for the correlation length is 



g a=^(l + ^) (17) 

and the function c(e) is of order one. The phase transition is at a = 1, which corresponds 
to Ej/Eq = 1/tt 2 . Thus, without vortex fluctuations the phase transition is at a smaller 
Ej/Ec value than the 2 /n 2 that follows from a duality argument [12]j. 
The correlation function is modified to 

(m)^ = Col = + ^E£k 2 le. (18) 

Thus the spinwave dispersion is affected at small distances (large k- vectors). This hardening 
of the spinwave dispersion close to the transition may be interpreted as resulting from the 
discreteness of the Cooper pairs, which makes fluctuations of phase and charge on short 
distances unfavourable. It leads to a mass 



M„ = — ^ In 



8tt? 



2f_e 

eE, 



1 + 



c 



(19) 



In the limit of small £ the Eckern-Schmid mass is recovered. An extrapolation to the S-I 
transition where £ — * oo yields a mass that vanishes at the transition, see Fig.2. We find a 
similar result for the instanton action 



Again, the classical result is recovered in the limit £ —>■ oo, whereas an extrapolation to the 
transition gives an instanton action that vanishes. From the WKB relation between mass, 
instanton action and barrier height Ub ar (as discussed in Sect.pTTD, we find that the depinning 
current Id ep ~ Sf nst /M v ~ l/ln(£) and thus vanishes algebraically close to the transition, 
i.e. Idep ~ (1 — a) 5 where a was defined in Eq. (|T7|) . 
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With the charge-charge correlation given in Eq. ([18|) we may calculate the spinwave damp- 
ing of vortex motion due to the coupling to spinwaves beyond the classical limit. Replacing 
ujk by Uk in Eq.(|12|), the overlap integral over the delta functions only contributes for vortex 
velocities that are higher than a threshold velocity vt = 2n^Ej/e (see Fig. 2). Note that 
this threshold velocity is independent of the momentum integration cut-off that is used. 
Thus, for vortex velocities v < v t there is no constant velocity solution to the equation of 
motion, unless the external driving current 1 = 0. Taking into account quantum effects 
changes the spinwave spectrum in such a way that the velocity window in which vortices 
move over the lattice potential without emitting spinwaves grows larger. The resulting re- 
lation between applied current and vortex velocity is shown in Fig.l for several values of vt- 
An extrapolation to the S-I transition yields a diverging threshold velocity and vortices and 
spinwaves decouple. 



V. SHORT RANGE COULOMB INTERACTIONS 

In this section we consider the case of short range Coulomb interaction U, i.e. the junction 
capacitance is negligible compared with the capacitance to the ground. This limit is of more 
relevance for 2-dimensional superconducting films. The system undergoes a T = phase 



transition which belongs to the 2+1-dimensional XY universality class [[□]]. The critical 
properties are well captured by a coarse-grained Ginzburg-Landau free energy. The effec- 
tive free energy has been derived from the Hamiltonian (§) using a Hubbard- Stratonovich 
transformation Jll , |27 |. 



It is therefore natural to express the charge-charge correlation function in terms of a 
Ginzburg-Landau coarse grained order parameter field. The charge-charge correlation can 
be expressed as a functional derivative of an appropriate generating functional as follows 



see Apppendix [B 1| for more details) 
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(qj T qkr') = j= — -= In 

OfIj T OfikT' 



T>if>T>if> exp(— F[if), if), //]) 



where 
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F=f<Px{^ |V^| 2 +r \ ^j\ 2 +u |^| 4 +C \ {d T -fi)^\ 2 } , (21) 

and the coefficients are r = l/2Ej - l/2E a , u = 7/128^, and C = 1/32^. The charge- 
charge correlation is related to the response of the system to a twist of the boundary condi- 
tions in the time direction. Its Fourier transform is 



2(V>(r, r)^*((r, r)) - 4( d 2 rdr(J T (f, r) J r (0, 0))e 



Where the current in the r direction J T is defined as 



ikr+iuiT 



(22) 



J r = i {^(f, r)<9W(r, r) - V(r, r)^*^, r)} (23) 

This type of correlation function was extensively investigated in [[TJJ. Due to the isotropy 
of the model in space-time the k = uj = term in Eq. (|22|) is proportional to the superfluid 
density p s of the system. 

As discussed previously the adiabatic mass is related to the zero frequency component 
of the charge-charge correlation function. Therefore one obtains 

M v ~ p s \n(L/a) + terms not divergent in L , (24) 

where the terms independent of L, that do not diverge with the system size, may arise from 
the /c-dependence of the zero frequency component of the charge-charge correlation function, 
however, the ln(L/a) is dominant in the superconducting case. Close to the transition 
p s ~ (Ej/E — with P « 2/3. We stress that this result is independent on the particular 
approximation we may choose to evaluate the correlation functions. 

More care it is needed right at the transition where it can be shown that the vortex 
mass does not vanish; we remind that it is determined by an integral of the charge-charge 
correlation function over the first Brillouin zone. Although (qq) k=LU=0 vanishes, there is an 



important contribution from the k dependence of the correlation function. Following [15] it 
may be calculated employing a 1/N expansion and to leading order 
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{qq)k ^ = 327E V k2 + UJ * /AEl (25) 
13 



The k-dependence at zero frequency will regularize the k-integral for the mass and, as a 
result it will not vanish at the transition but becomes independent of the system size 

7T 3 / 2 

M trans = — — « 0.44M o (26) 

In the superconducting phase the charge-charge correlation function can be approximated 
to 

{qq) k , u „~p s k 2 /{k 2 + A(ujl) (27) 
Spinwave damping may be calculated from the equation of motion. From Eqs.(|l|) and (27) a 



threshold velocity v t = Cu p y E /8Ej is found, which for the short range Coulomb interacting 
does not diverge at the transition. 

We checked the main conclusions of this section performing simulations that are presented 
in the next section. 



VI. MONTE CARLO SIMULATIONS 

We now turn to the Monte Carlo results that allow for a check of the previous calculations 
and provide information for the region where the self consistent harmonic approximation for 
the long range case is not valid. Since the CCVG described by Eq. @ contains an imaginary 
coupling and long range interactions, it is more convenient [[U| to simulate the system in the 
equivalent current loop representation Eq . (|A5|) . The correspondence is simply (qq) = (J°J°). 
The condition that the currents J M be divergenceless is taken into account by making Monte 
Carlo steps that preserve the property V^J^ = 0. Thus, we create or annihilate small 
current loops, as well as 'periodic current loops' that go through the whole system which is 
taken to have periodic boundary conditions in all three directions. In the case of long range 
Coulomb interaction periodic current loops in the time direction are forbidden since they 
violate charge-neutrality as demanded by the logarithmic Coulomb interaction. 

The size in the time direction was taken to be 8, which corresponds to a temperature T 
equal to one eight of the plasma frequency, which is quite low for a JJA. The simulations 
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were done on lattices of linear dimension varying from 4 to 12 and the standard Metropolis 
algorithm was used, with typically 5000 sweeps through the lattice for equilibration and the 
same amount for measurement. In all cases the full charge-charge correlation function was 
measured and with the help of Eq.(|8|) the mass was determined. Reliable results for the 
correlation length would require larger systems than we were able to simulate. 

For logarithmic Coulomb interaction data for the vortex mass is shown in Fig. 3. The 
instanton action Si ns t for a vortex hop from a plaquette to a neighboring one and therefore 
the depinning current behave in a similar way. The critical point is at Ej/Ec ~ .6, which 
agrees with experimental findings . Close to the S-I transition the vortex mass depends 
strongly on the ratio Ej/Eq- Note that apart from small corrections the mass is system size 
independent. 

The results for the vortex mass for onsite Coulomb interaction are shown in Fig. 4. The 
critical point is at Ej/E Q ^ .85. The collapse of the curves for different system sizes (see 
Fig. 4b) demonstrates that the vortex mass indeed scales with ln(L) in the SC phase, whereas 
in the insulating phase it scales to zero (see the inset of Fig.4b). If the logarithm of the system 
size is not scaled out, the curves for the mass in systems of different size approximately 
cross at the transition between 0.4 and 0.5 times M , which is in good agreement with 
M trans /M m 0.44. The region were the mass is strongly dependent on the ratio of the 
couplings is somewhat larger than in the long range case. 

VII. CONCLUSION 

We have investigated vortex motion, mass and spinwave damping in Josephson junc- 
tion arrays. We first derived a one-vortex effective action from which it became clear that 
dynamical vortex properties are governed by the charge-charge correlation function. We 
showed how to recover all known results for classical arrays. The Eckern-Schmid mass being 
proportional to the junction capacitance, was found to be correct in the limit Ej 3> Ec, but 
also the ground capacitance contributes to the classical mass, which becomes system-size 
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dependent. 

In the generic situation vortex motion is affected also by the presence of other vortices 
(or more precise dipoles in the superconducting state) which are present because of thermal 
or quantum (this is the case we considered in detail) fluctuations. In the quantum regime 
close to the S-I transition we investigated the vortex properties analytically and by means 
of Monte Carlo simulations. 

We considered specifically the two extreme cases in which either the self-capacitance C D or 
the junction capaciatnce C was set two zero. In the case of long range Coulomb interactions 
(C equal to zero) the main conclusions are that the mass and depinning current vanish at the 
phase transition in a way that reflects the nature of the S-I transition, whereas the velocity 
window in which vortices can move without exciting spinwaves grows. Our predictions for 



the depinning current are in qualitative agreement with experiment [ 29| . Our results suggest 
that ballistic vortex motion may be seen best in arrays that are close to the S-I transition. 

The limit of on-site Coulomb interaction seems more appropriate for granular or uniform 
films. In this case the mass has a logarithmic dependence on the size of the sample in the 
superconducting region. It is proportional to the superfluid density, and therefore shows 
critical behaviour approaching the superconductor-insulator transition. At the transition it 
does not vanish, but becomes scale independent. 

Using Monte Carlo simulations for determining the charge-charge correlation function 
numerically, we were able to verify several conclusions of the self consistent harmonic ap- 
proximation (for long range Coulomb interactions) and the coarse-graining approach (for 
short range Coulomb interactions). 

In all the results we presented we were mainly concerned with the superconducting side 
of the S-I transition, allthough the equations derived in section are in principle valid 
throughout the phase diagram. We might also investigate the dynamical properties of the 
vortices in the insulating region. This is, however, not useful since the vortices are massless, 
delocalized, and strongly fluctuating in the insulating phase. In particular when E c ^> Ej 
the charge on the superconducting islands is a good quantum number and it is more useful to 
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describe the quantum dynamics of charges in the resistive (or insulating) region, which may 
be investigated using the same techniques as presented here for vortices. In the resistive 
(high temperature) phase the concept of vortex is also not usefull as the system is not 
globally superconducting and the superconducting phases of the islands are disordered. We 
may also add that the main approximation that leads to the vortex effective action, namely 
the truncation of the cumulant expansion to the terms quadratic in the velocities may not 
be justified in the disordered (both resitive and insulating) phases. Probably in this case a 
more complicated equation of motion that includes terms proportional to higher powers of 
the velocity should be included. 
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APPENDIX A: DUALITY TRANSFORMATIONS 

This appendix reviews shortly how to pass from Eq.(0) to Eq. ([$]). For more details and 
the treatment of the external current contribution we refer to Ref. |TJ|. We start from the 
basic expression for the partition function Z 

Z = Tr exp(-/3#), (Al) 

where (3 is the inverse temperature T and H the Hamiltonian Eq.([2|). We go over to a 
Euclidean path-integral formulation by introducing time-slices, i.e. dividing /3 in N T intervals 
of size e, such that N T e = (3. Inserting complete sets of states at each time slice we arrive at 
f I 1aE ■ 

Z = J Vfar exp < qirUijqjr + i J2 QiT&T + eE J Yl COS (^r ~ 4>jr) \ (A2) 

{q iT } { ^ ijr ir <ij>r ) 

In Fourier components, the electrostatic interaction between the charges is Uk = 27r/ (/c 2 +A 2 ), 
where the inverse range of the interaction is A = yj C j C. It is related to the inverse 
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capacitance matrix by Uij = 2-nCC^ 1 . Now we make the Villain approximation [30] for the 



cosine term 



^-f(eEj)(cf) iT - <j) jT - 2ixn ij)T y \ , (A3) 



eEj 

exp i - 

where n is a directed discrete field that lives on the bonds between lattice sites. The function 
/ equals unity if its argument is large, i.e. when the product eEj is not too small. After a 
subsequent Poisson resummation, i.e. writing 

E F[n ijjT ]= Y J T>nF[n ijjT ]exp(2Tri^2n ij:T J ijiT ), (A4) 

an integration over the fields n%j 7T and the phases <fri T yields a representation in terms of 
divergenceless discrete current loops 



Z= E <W) exp £ JlU tJ Jl - — J2(Jtr) 2 , (A5) 

{J? T } { 71 H,r ^Ej ira 



where a = x,y and \i — x, y, 0. Here the time components J° of the current are simply 
the charges. The vortex degrees of freedom may now be extracted by solving the constraint 
in Eq . (|A5|) by writing Jf T = e ab VbipiT — O a qi T (the operator O a denotes the line integral in 
direction a, i.e. it is the inverse of V a ) and making a final Poisson resummation on the 
discrete field ip. Making use of the identity 0^ = e ab V a O b Gij Eq.(|3p results. 

Note that we did not keep track of determinants, as they are irrelevant for the present 
purpose. The requirements that both the factorization of exp(—(3H) and the Villain ap- 
proximation are valid restricts e to be of the order of the inverse of the plasma frequency 



uj p = y/8EjEc for long and uj p = \/%EjE for short range Coulomb interactions, i.e. we 
take e = p/oo p or e = p/oo p . Details and numerical factors may depend on the exact choice 
for p. We take p — 1, except for the self consistent harmonic approximation where we take 
p = v^2- The plasma frequency is the natural frequency for spinwaves. 
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APPENDIX B: THE CHARGE-CHARGE CORRELATION FUNCTION 



In general the charge-charge correlation function may be expressed as a functional deriva- 
tive of the free energy in the following way 

5 2 



(Qjtqw) 



In 



(Bl) 



At=0 



Here £ denotes the field or fields that are integrated (or summed) over, S[£] a corresponding 
action and / the appropriate measure. 

1. long range Coulomb interactions 

For instance, if one takes the action to be the charge part of the CCVG action, as is 



done in section |V|, one has 

8 2 



(qktqw) 



In 



if 



X] eX P - qirQijrr'qjT' + 



:b2) 



Qir \ IJTT' 

The kernel Q was defined in Eq.(^). Now the partition function may be rewritten in terms 
of different (dual) fields and if we keep track of the 'currents' \x during the transformation 
we may express the charge-charge correlation function in terms of the correlation function 
of the new fields. Applying this strategy to Eq.([B^), we find after a Poisson resummation 

5 2 

— 111 I , 



(qktqif) 



Sfi kt Sfx, 



-In [j>Xp ( - 7T 2 UrQi^Jjr' + 



ijrr' 



+Z7T J2 hrQijlr'^jr' + ^rQ lj l T ^j r 



ijrr r 



ijrr' 



(B3) 



from which we read of immediately that 



Iktqif) — -jQkltt' ~ 71-2 Y QkrntT^mrlnT^Qnlrit'- 



(B4) 



This has simplified the problem of calculating the correlation function considerably, since 
the new field I interacts with the kernel Q -1 which is, in contrast to the original kernel Q, 
short range in both the space and time directions. This discrete Gaussian model is therefore 
convenient for Monte Carlo simulations. 
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Another representation for the correlation function is found by applying an inverse Villain 



approximation, see chapter 11 of Ref. |25[ for details. This enables one to rewrite the 
partition function as a path integral over a continuous field with sine-Gordon action. The 
result is very similar to Eq. 



(qjrqkT*) = ^Qjkrr' ~ "^^Qjm 



mrtWintiPnt' 



)Qnkt>T'i 



(B5) 



mnW 



but now the correlation function of the dual variables is now to be calculated using the 



action 25 



Sty] =7i 2 J2Yl fctQijtffa* ~ H Y1 cos(2vr0,i) 

ij tt' it 



(B6) 



In this formulation we may employ the self consistent harmonic approximation |26| 



2. short range Coulomb interactions 



Here we outline the derivation of Eq . fl2~T|) in the text. For details we refer to Ref. f27 |. In 
the case of short range interaction, the critical properties of the system are well described 
by a coarse grained free energy. The generating functional that we consider is the parti- 
tion function Eq.flX^) with a coupling to a current \x as in Eq. (|B"T|) . After performing the 
integration over Dq we have in terms of phase-variables only 
8 2 



(QjtQkt') 



SfJ-jtSfJ-kt' 



In 



ij 



~Ej E / dT 



'COS 0,; - 



, (B7) 



where the measure T><fi still contains a summation over winding numbers (i.e. 4>i(/3) = 
0i(O) + 27rmj) in order to account for the discreteness of the charges. Decoupling the Joseph- 
son term by means of a Hubbard- Stratonovich transformation one gets: 

5 2 . - / I I 



(QjtQkt') 



8ii jt 8ii 



In 



kf 



VipVip exp 



2E, 



eX P i E / dr \^i( T ) e 



[W) 
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where the matrix t < ij > is one for nearest neighbors and zero otherwise. The average in the 
last factor is with respect to the remainder of the action in ( |B7| ). Close to the transition 
point we may expand the average in powers (cumulants) of the fields ip; this expansion yields 
the generating functional used in the paper. We note that the currents \x enter the time 
derivatives in the generating functional Eq.(^TJ) in a gauge invariant way. 
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Figure Captions: 

Fig. 1 : Driving-current vs. vortex velocity relations for long range Coulomb interactions 
and several couplings. From left to right: v t = (for Ej/Ec 3> 1, the classical limit, see 
0), Vt/uJp = 0.5 and Vt/u p = 1. The vertical dotted line indicates the minimum velocity a 
vortex needs to move ballistically over the pinning potential. 

Fig. 2 : Mass, threshold velocity and coherence length as from the self consistent harmonic 
approximation. 

Fig. 3 : Monte Carlo result for the mass with long range interactions. Shown are a 6x6x8 
(crosses) and a 10x10x8 (diamonds) system. The statistical errors are of the order of the 
symbol size. 

Fig. 4 : Monte Carlo results for the mass with on-site interactions, a) Shown are a 4x4x8 
(lower curve), 6x6x8, 8x8x8, 10x10x8 and a 12x12x8 (upper curve) system. The statistical 
errors are of the order of the symbol size. 

b) The mass for on-site interactions with the logarithmic of the system size scaled out. 
Clearly the curves collapse on one line in the SC phase, but not in the critical region. Inset: 
The size dependence of the mass with the logarithm scaled out for Ej/E Q =1.45, 1.1, .95, 
.85, .75 (from top to bottom). 
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